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We study the lepton flavor violating processes of r^ — > n^PP decays with PP — 

K'^K~ , K'^K'^,tv'^tt~ ,Tv^TT^^ in the framework of two Higgs doublet model III by virtue of the Chi- 

^—^ ral Perturbation Theory. In this model, only three neutral Higgs bosons contribute to these de- 

IT^ cays. With the current experimental constraints, we show that (a) the contributions of \\uu{dd)\ 






X 



term are very small for these four decays; (b) we get the correlation between \Xas\ and |Arp|. 
For \\r)i\ ~ 10 — 400, one has \\ss\ ^ 40—1; (c) in the existing model parameter space, 
Br(T^ — > ji~ K^K~) could reach the order of C>(10~*), but Br{T — )> ij,~-k^-k~ /ti'^-k'^) are too 
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^^ Flavor physics has made rapid development in these decades. In addition to B physics, r physics, including 

^^ the determination of as from the inclusive hadronic width, tests of charged-current universality, and lepton flavor 

O 

'""^ violation (LFV) decays, etc., also belong to one branch of flavor physics. Among these topics, LFV decays became 

> 

more important after the discovery of neutrino flavor oscillations and related non-zero neutrino masses [TH3] . The 



^ LFV decays are forbidden or suppressed strongly in the Standard Model (SM). While for new physics, their decay 

ratios would be enhanced largely by new LFV sources and/or new particles effects. Obviously, LFV decays would 
be a sharp tool to seek for the new physics beyond the SM. Due to its large mass, r lepton has many LFV decay 
modes, such as r — >■ 31', r — > I'j and r ^- l'M{M — P,V^ ,PP), which could display rich information on LFV 
mechanism. In particular, in certain new physics models, e.g., where the LFV effects arise through the Higgs boson 
exchanges, the t lepton is expected to have large LFV couplings to the /i because of its large mass. Therefore, t 
LFV decays become an interesting flcld for finding new physics. 

Among r various LFV decays, r — ?► 2>l' ,t — ?► n^y and r -^ I' PiV'^) have been searched in experiment j4, and 
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investigated in many theoretical frameworks, such as reahstic non-minimal supergravity models [5], SUSY SO(IO) 
models [3 E] , various extended-MSSM models [HI [S] , Seesaw models [TUHIl] , Unparticle scenario [T3] , Topcolor 
models [H], the Littlest Higgs model with T-Parity [iS], etc. For t^ -> fjTPP decays, the newest experimental 
upper limits of PP = K+K-{K"K", tt+tt- ,n^7r") decays are [1[I1[I7]: 

B{t- -^ H-K+K-) < 6.8 X 10~^ 90%CL , 
B{t' -^ /i-TT+TT-) < 3.3 X 10-^ 90%CL , 
B{t- -^ fi-K°K"^) < 8.0 X 10"^ 90%CL , 
B(t- -^ /x-7r"7r°) < 1.4 x 10~^ 90%CL , 

where the first three values have been improved by roughly one or two orders of magnitudes. Accordingly, their 
theoretical investigations have been done in many possible extensions of the SM |jT5ti22j . C. Chen et al. have 
analyzed the scalar bosons effects on t^ — > /.i^ K^ K^ {K'~' K'^) decays in the supersymmetric seesaw model with 
nonholomorphic terms in the lepton sector at the condition of large tan /3 [19] . E. Arganda et al. have investigated 
these processes in two constrained seesaw scenarios in the Minimal Supersymmetric Standard Model (MSSM) [20] . 
M. Herrero et al. have found that t^ -^ ji^-q and r^ -^ A* /o are more efficient to test indirectly to the Higgs sector 
than T — > 3/i,T^ — > ji^K^K^ decays in supersymmetric seesaw models [21]. And C. Yue et al. has calculated the 
new particle effects in the topcolor-assisted technicolor model and the littlest Higgs model with T parity I22| . 

In this paper, we shall study r decays in the Two Higgs Doublet Model HI (2HDM HI) since it naturally introduces 
flavor-changing neutral currents (FCNCs) at tree level. In order to satisfy the current experiment constraints, the 
tree-level FCNCs are suppressed in low-energy experiments for the first two generation fcrmions. While processes 
concerning with the third generation fermions would be larger. Therefore, these FCNCs with neutral Higgs bosons 
mediated may produce sizable effects on the t — fi transition. Moreover, the evaluation of neutral Higgs effects 
would help to provide useful information to identify the Higgs signals in LHC experiment. Hitherto, many works 
have been performed in this scenario. For example, the authors in [23j have discussed the r — > 3/i decay in the 
2HDM HI and MSSM, and analyzed their collider signals, r* — > Z^P" decays have been considered in the 2HDM 
HI model with four-tcxturc Yukawa couplings [24\ , which is shown that the model parameters in leptonic sector are 
less constrained by present experimental data, r — > ^P{V°) decays have been studied in our previous paper [^[^ . 
where the hadronization in final state is merely expressed in terms of the meson decay constants and meson masses. 
Furthermore, t ^ fi transition mediated via Higgs bosons is sensitive to the Higgs sector. In theory, the Higgs 
influences on LFV decays have argued in 2HDM HI and MSSM [?7H5n] . Aoki et al. 31 have disputed the decay 
properties of Higgs bosons in four types of 2HDMs and their collier phenomenology. So we extend our discussion 
to the case of two pseudoscalar mesons in the hadronic final state and deal with hadron matrix elements by the 



chiral perturbative Theory (^PT) 

The paper is organized as follows. In Section II, we make a brief introduction of the theoretical framework for 
the 2IIDM III. In Section III, we deliberate the calculations of the decay amplitudes with xPT and our numerical 
predictions. Our conclusions are given in the last section. 

II. THE TWO-HIGGS -DOUBLET MODEL III 

As the simplest extension of the SM, the Two-Higgs-Doublet Model has an additional Higgs doublet. In order to 
ensure the forbidden FCNCs at tree level, it requires that either the same doublet couple to the w-type and d-type 
quarks (2IIDM I) or one scalar doublet couple to the u-type quarks and the other to d-type quarks (2IIDM II). 
While in the 2HDM III [SSHIO] , the two Higgs doublets could couple to the w-type and rf-type quarks simultaneously. 
Particularly, without an ad hoc discrete symmetry exerted, this model does have FCNCs at the tree level. 

The Yukawa Lagrangian is generally expressed as the following form 

where Hi{i ~ 1,2) are the two Higgs doublets. Qi^L is the left-handed fermion doublet, Uj,u and Dj^r are the 
right-handed singlets, respectively. These Qi,L, Uj_r and Djji are weak eigenstates, which can be rotated into mass 
eigenstates. While rj^'^ and £,^'^ are the non-diagonal matrices of the Yukawa couplings. 
We can conveniently choose a suitable basis to denote Hi and H2 as 

\ f V2G 

v + ^l) \ tG" 

where G^'^ are the Goldstone bosons, H^ and A" are the physical charged-Higgs boson and CP-odd neutral Higgs 
boson, respectively. Its virtue is that the first doublet Hi corresponds to the scalar doublet of the SM while the 
new Higgs fields arise from the second doublet H2. 

The CP-even neutral Higgs boson mass eigenstates H^ and /i° are linear combinations of (/)" and (j)2 in Eq. ^ 
as follows 

H'^ = 0° cosa + (/)2 sina, h'^ = —(f)lsma + (jy^cosa, (3) 

where a is the mixing angle. 

After diagonalizing the mass matrix of the fcrnrion fields, the Yukawa Lagrangian becomes |41j 
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where U and D are now the fermion mass eigenstates and 
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where Vj^ '^ are the rotation matrices acting on up- and down-type quarks, with left- and right-chirahties respec- 
tively. Thus, VcKM = {yL)^^L is ttis usual Cabibbo-Kobayashi-Maskawa (CKM) matrix. In general, the matrices 
YjU,D jj^ gq |g| g^j,g diagonal, while the matrices S,^'^ are non-diagonal which could induce scalar-mediated FCNC. 
From Eq. Q, we obtain that the couplings of neutral Higgs bosons to the fermions could generate the FCNCs. 
For the arbitrariness of definition for ^- ■' couplings, we can adopt the rotated couplings expressed ^^'^ in stead 
of ^'^•^ hereafter. 

In this paper, we use the Cheng-Sher ansatz |35j 



t-U.D , v2y^ 
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which ensures that the FCNCs within the first two generations are naturally suppressed by small fermions masses. 
This ansatz suggests that LFV couplings involving the electron are suppressed, while LFV transitions involving 
muon and tau are much less suppressed and may lead to some effects which are promising to be tested by the 
future B factory experiments. In Eq. (|7|, the parameters \ij are complex where i, j are the generation indexes. In 
this study, we will discuss t^ — > fi^PP decays in the 2IIDM III. 

III. THE DECAY AMPLITUDES FOR r -^ /I-PP DECAYS 



In 2IIDM model III, due to the existence of FCNCs, the t^ -^ ji^ PP processes could occour at tree level with 
the neutral Higgs bosons H^^hP and A° mediated. The decay amplitudes could be factorized into leptonic vertex 
corrections and hadronic parts described with hadronic matrix elements, which be expressed as: 



in PP\M\t ) = —j^-mq^/m~w^- 



H'^* ■ A.^ • (Ms+p + H'^ ■ Kn ■ (Ms-p 



<PP\iqq)s\0> 



N''* ■ A.^ • {flT)s+P - N'i ■ A* • {nT)s-i 



<PP|(q-q)p|0> , 



(8) 



where mq{q = u, d, s) are the masses of quarks which constitute the meson in the final states, mr and tti^, are the 
masses of lepton t and fi. S ± P denotes (1 ± 75). The auxihary functions H'^ and N'^ can be written as 

ff" = (A + B)- A:„ + iA-B)- \uu + F, H''-^ ^{A-B)- A*,(,,) + {A + B)- X^iss) + F, 

iV« = {A + B)- A:„ -{A-B)- A„„, TV^'^ ^{A-B)- A:,(,,) -(A + B)- A,,(,,), 

sin^ a cos^ o; ^^ 1 ^^ r, • / 1 1 \ /«\ 

A = — ^ 1 5 — , B = — ^— , t ~ 2smacosa( — 7, s-j, (9) 



where mno, m/jo and Wyio are the masses of neutral Higgs bosons, a is the mixing angle. A^-^t is the LFV coupling 
parameter of neutral Higgs bosons to lepton r and /i, and Xuu{dd.ss) are the quark counterparts. < PP\{qq)s{p)\0 > 
are the hadronic matrix elements. From the Eq. (Is]) and (lo]), we can see that three neutral Higgs bosons take effects 
to two pseudo-scalar mesons through {qq)s±p operators. But it differs from the MSSM where the 7 contributions 
are the dominant one and only H'~' and h'^ take effects at the case of large tan /3 [5D] . 

Next, we will calculate the hadronic matrix elements. The ordinary method is to parameterize the hadronic 
vertexes by hadron masses and their decay constants. For the r^ — > fi^PP decays, its scale lies in the non- 
perturbative region and vector and scalar resonances generally play a propelling role. However, the masses of 
scalar resonances are so high that their effects are ignored. Besides, the heavy Higgs bosons cause the hadronic 
final state insensitive to resonances and their influences are less known. Therefore, we apply xPF [32H34] instead 
of Resonances Chiral Theory {RxT) [HI |12I to handle < PP\{qq)s{p)\Q >■ Additionally, it should be noted that 
the pseudo-scalar currents have contributions only to one pseudo-scalar meson in final states. Thus, by virtue of 
the relevant scalar currents S'*(i — 0, 3, 8) in xPF and using the following currents [5D]: 

-uu = ^^3+ 1 58^ 1 50 
2 2V3 VG 

-dd = -^S^ + ^S' + ^S°, 
2 2^/3 VQ 

-ss = -^S^ + ^S°, (10) 

we could obtain the amplitudes: 

M{T-^fI-pP) = '^■^m^^-[T{PP)-Xr^-{fiT)s+P+T*{PP)-X;^-{flT)s-p], (11) 

T{K+K-) = A ■ [Re{Xuu) ■ ml + Re{X,,) ■ {2ml - ^1)] + Bi ■ [/m(A,,) • (2m| - ml) - /m(A„„) • ml] 
+F-ml, (12) 

r(if"X°) = A ■ [Re{XM) ■ ml + Re{X,s) ■ {2ml - ml)] + Bi ■ [/to(A,,) • {2ml - ml) + Im{Xdd) ■ ml] 

+F-ml, (13) 

T{iT+n-) = ml ■ {A[Re{Xuu) + Re{XM)] + Bi ■ [Im{XM) - /m(A„„)] + F} , (14) 

2 

n^V) = ^■{A[Re{X^^)-Re{Xad)]-Bi-[Im{Xuu)+Im{Xdd)]}, (15) 



where TOtt and rriK are the masses of tt meson and K meson. The expressions of A, B, and F can be found in 
Eq. (|9|. The T{K^K^) and T{K'^K^) terms are relevant to model parameters \uu{ss) and Xdd{ss): respectively. 



The differences of Eq. ( 13 1 from Eq. ( 12 ) are that Xdd replaces A„„ and there are an opposite sign before Im{\dd)-m^ 



term. While both the T(7r+7r ) and r(7r°7r°) terms only get involved in the parameters A„„ and A, 



dd- 



IV. NUMERICAL RESULTS 

In our calculation, the input parameters are the Higgs masses, mixing angle a, and couplings \qq{q = u, d, s), A^-p, 
and their phase 9. We use the values of parameters in the literatures [H] [IS] 



niHO = leOGeF, m^o = 115GeF, m^o = l2QGeV, a = 7r/4, 6*,, = 7r/4, 



(16) 



where the experimental constraints from B — B mixing, b — > S7, p^ , and i?;, are considered. The constraints 
on \\uu(dd,ss)\ and |At-^| from different phenomenological considerations |45H51j are demonstrated in Table I. The 
constraints on \Xss\ from Bg — B^ mixing and h ^ s^ process, are ^0(1) ;45J. Whereas for the first generation, 
the FCNC couplings are suppressed, the values of Xuu and Xdd are less than 1 [55]. For parameter in lepton sector, 
the r -> 3^, T -> ^7, h° — > // decays and {g — 2)^ have constrainted A^^ at the order of 0(10) — O(IO^) jiMST] . 
Hence, for the following values^*! of \Xrfj\ = 150, |A,(„(dd)| ^ (0, 0.9), \Xss\ ^ (1, 5), the decay ratios are calculated. 
First, our calculations show that these branching ratios are dependent on ^t-^, the phase angle of Xr^- As in 
literatures, we assume the value of 6',-^ as tt/A. 

TABLE I: Constraints on the A^j in quark and lepton sectors. 





Bounds and restrictions 


Process and Restriction 


References 


Quark sector 


|As.| >C»(1) 


Bs — Bs mixing, 6 — >■ 37 


m 




\Xuu\, Add ~ 0{l) 


F" - FO mixing (F = K, Ba, D), R^, p,B ^ X,y 


m 




A,„ ~ 0(10) - 0(10") 


(g-2)^,mA0 — > 00 


[49]** 


Lepton sector 


A^^ ^ Arr ^ A^T ^ Ae/j ^ 10 


h"' ^ ff 


EQ] 




A,^ ~ C»(10") - C»(10') 


T — > 3/1, r — > /i7 


m** 



Note that the constraints in i49i and |5T] are denoted by our notation. 



The Fig. 1 shows when \Xt-^\ ~ 150, the functions of Br{T -^ p KK) versus model parameters, where (a) is for 
T^ — > fj.^ K^K^ decay and (b) is for r" — > p^ K'^K^ decay. One can see from Fig. 1(a) that Br{T^ — > ijT K^ K^) 



[*] The constraint on |At^| is found to bo |At-^| < 2.76 in Ref. '25 ,. which is because there the values of \\uu\ = 150 and |Add| = 120, 
completely different from those used in this paper, are used. 



raises with the increase of \Xss\ but does not vary with |A„ii| growing. For t~ -^ ^i~K^K^ decay, the same is as 
that of T^ — > fi^K^K^ decay except that \Xdd\ replaces |A„,i|. In the mentioned parameter space, their decay 
ratios could reach the experimental upper limits. Thus, with fixed value of |At-^| , Br{T^ — > ^T KK) are insensitive 
to the variation of lAuurrft;-)! and the contributions of jA^sl are larger than those of |Au„((j£i)|. As can directly be seen 
from the Eq. (12) and (13), these results are expected because the squared mass of K is significantly larger than 
the squared mass of pion unless the value of |A„„(rf£()| is significantly larger than \\ss\- 

With the upper bounds of Br(T^ — > ii^ KK), we take |A„ii| — \Xdd\ = 0.5, and draw the constraint on lAr^^l and 
|Ass|. Their correlation is shown in Fig. 2, where the dependences of lA^-^l on \\ss\ are denoted by the solid line 
for r^ -^ yT K^K^ decay and the dotted line for t^ -^ [iT K^iC^ decay, respectively. One could get from Fig. 2 
that, under the current experimental bounds, when |Ar;^| becomes large, the \\ss\ goes to decrease, and vice versa. 
For example, the LFV coupling |At-^| — 10, the upper limits of Br(T^ — > fi^K^K^) and Br{T^ — > n^ K^K'^) on 
\Xss\ are at 40 and 60, respectively. When |Ar^| = 400, the values of |Ass| are bounded at 1 for t^ -^ fi^K^K^ 
decay and 1.4 for t^ -^ ji^K'^K'^ decay. Obviously, the upper bound from Br{T^ -^ fi^K^K^) on \Xss\ is more 
stringent than that from Br(T^ -^ fi^K'^K'^). And in the area under the solid curve, all the values of \Xr)j,\ and 
I Ass I are allowed. 

t-->mA'+r t-->mK'>K° 
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FIG. 1: For jA^pl = 150, the functions of branching ratios versus model parameters. Here, (a) is for r^ — >■ y" K'^ K~ decay, 
while (b) is for r^ — > y~ K^K'^ decay. 

The functions of Br{T^ — > ii^tttt) versus |A,i„| and \Xdd\ are presented in Fig. 3, where [A,-^! = 150, (a) is for 
T^ -^ /_t^7r+7r~ decay, and (b) is for r^ — > /i^Tr^vr*^ decay. In Fig. 3, with IAt-^I = 150, Br{T^ — > /i^Tr+Tr^/vr^Tr") 
move upward with the rising of |A„„| and \Xdd\, which is different from the case of r^ — )■ fi^ KK decays. But their 
decay ratios only reach to the order of 0(10^^"'^). Even when |At-^| — 450, the values of Br{T^ — > /i^Tr+vr^/Tr^Tr") 
are only improved one magnitude of order. Apparently, in order to reach the experimental data, |At-^| should be 
enhanced to 0(10'^) and more. Furthermore, comparing to Fig. 3(b), the branching ratio of t^ — > /i^7r+7r^ decay 
keeps a relatively rapid growth and is more sensitive to |A„„| and \Xdd\- 
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FIG. 2: The constraint on |At-^| and \Xss\ with the experimental upper bounds on Br{T — >■ fi K'^K ) and Br{T 
li-K°K") which read Br{r- -> ^.-K+K') < 6.8 x lO'^iCL = 90%CL) and Br{T' -^ ^- K° K'^) < 1.6 x 10"^(CL 
90%CL). 
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FIG. 3: For |At,i| = 150, the functions of branching ratios versus model parameters. Here, (a) is for r^ — > h^tv^tt^ decay, 
while (b) is for r^ — > /^^tt^tt" decay. 

V. CONCLUSION 



In this paper, we have calculated the branching ratios of r^ — > ^i^ PP{PP = K^K^ ,K^K'^, tt+tt^, 77*^77") decays 
in the 2HDM III. Because there are FCNCs at the tree level, all three neutral Higgs bosons have mediated theses 
decays by scalar and pseudoscalar currents. This is different from the MSSM where only H'^ and h'^ play roles in 
these processes. Since pseudoscalar currents have no contributions to two pseudoscalar mesons, we only considered 
the case of scalar currents. Although resonances play a role in t^ -^ fi^PP processes, the massive Higgs bosons 
are not sensitive to them. Consequently, we disregarded the vector resonance and used xP'^ to handle the hadron 
matrix elements. Our results suggested that, the |Au„(£;(i)| terms contribute really tiny to Br{T^ -^ fi^PP) due 
to the constraints on |Ai,„(dd)| from known experiments. As expected, Br{T^ — > ^i^ K^ K^) could give a stringent 
constraint on the parameters \\ss\ and \Xt^i\- We have used its experimental upper limit to obtain the correlation 



between \Xss\ and |At-^| and consequently the corresponding allowed region of \Xss\ for the given region of |Arp|- 
In the suitable parameter space satisfying all known constraints, Brir^ -^ /i^K^K^) could be up to 0(10^^). 
While for r^ — ^ fi^mr decays, the branching ratios are too small to be observed. It should be noted that the latest 
experimental values are based on 671 fb~^ data. We expect that the run of Super-B factory provides more and 
more data of LFV r decays and consequently to obtain the more stringent constraints on the model parameters. 
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